Surface vs. bulk Coulomb correlations in photoemission spectra of perovskites 
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Recent photoemission spectra of the perovskite series Sr^Cai-^VOs revealed strong modifications 
associated with surface contributions. To study the effect of Coulomb correlations in the bulk and 
at the surface the quasi-particle spectra are evaluated using the dynamical mean field theory. It is 
shown that as a result of the reduced coordination number of surface atoms correlation effects are 
stronger at the surface than in the bulk, in agreement with experiment. 
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Photoemission is a key spectroscopy for the investiga- 
tion of electronic properties of strongly correlated ma- 
terials. Because of its surface sensitivity there is grow- 
ing concern to what extent the data are representative 
of bulk properties and whether they need to be corrected 
for surface effects influencing the one- and many-electron 
properties. To study these questions Maiti et al. B and 
Sekiyama et al. || performed photoemission measure- 
ments on Sr^Cai^xVOs using a wide range of photon 
energies. As a result of the frequency dependent escape 
depth of the photoelectron, the spectra reveal striking 
variations evidently associated with surface modifications 
of the electronic structure: emission from the coherent 
peak near the Fermi level is reduced and correlation- 
induced satellites have larger weight than in the bulk. 
The same trend is observed on SrRuC>3 ||. Another ex- 
ample is S^RuCU for which earlier photoemission data 
appeared to be in conflict with de Haas- van Alphen mea- 
surements . Recent experimental and theoretical work 
proved that this discrepancy can be resolved by account- 
ing for the surface reconstruction of Sr2Ru04 which leads 
to significant changes in photoemission spectra || . 

The understanding of surface effects in photoemission 
spectra of transition metal oxides is important in order 
to distinguish them from correlation phenomena in the 
bulk. Although the electronic properties of strongly cor- 
related materials are currently a field of intense experi- 
mental and theoretical investigation || jjj, |) surface ef- 
fects have so far received little attention M . For this pur- 
pose photoemission studies of the series Sr^Cai-zVOs 
jj], ||, [DJ [l2| are very useful since the one-electron 
properties are relatively simple and the ratio U/W can 
be tuned systematically. Here, U is the on-site Coulomb 
energy and W the width of the relevant bands near the 
Fermi level. While SrV03 is considered to be a mod- 
erately correlated metal with a clear quasi-particle peak 
near Ep, the slightly smaller band width of CaVC-3 is be- 
lieved to push this system closer to the metal-insulator 
transition. For this reason CaVC>3 has become a key sys- 
tem of investigation |l3|, [l4|, For instance, recent de 
Haas-van Alphen measurements on CaV03 p5| detected 
a cubic Fermi surface although this system exhibits or- 
thorombic distortions. 

In the present work we use the dynamical mean field 



theory (DMFT) @, ||, to evaluate the surface and 
bulk quasi-particle spectra of SrVC-3. We show that be- 
cause of the planar electronic structure of perovskite sys- 
tems and the concomitant narrowing of the surface lo- 
cal density of states correlation effects at the surface are 
more pronounced than in the bulk. In the perovskites 
this enhancement is particularly strong since the local 
Coulomb energy is not far from the critical value for a 
metal-insulator transition. Such a trend had been pre- 
dicted by Potthoff and Nolting [§ who studied the metal- 
insulator phase diagram for a semi-infinite simple cubic 
s band at half filling. Here we calculate the self-energy 
for a multi-band system using a realistic local density of 
states appropriate for SrV03 and find qualitative agree- 
ment with photoemission data [[l], 

Electronic structure calculations for SrV03 within the 
local density approximation (LDA) |Q show that the 
conduction bands near Ep consist of degenerate ti g 
bands derived from V 4+ (3d 1 ) ions. The filled O 2p bands 
are separated from the ti g levels by a gap of about 1 eV, 
and the cubic crystal field of the V-0 octahedron shifts 
the e g bands above the ti g bands. The ti g bands can be 
represented via a tight-binding Hamiltonian with diago- 
nal elements h X y iX y(k) = erf, + to(C x + Cy) +t\C X Cy + \t% + 

h{c x + c y ) + t 4 c x c y ]c z , where a = 2 cos(fc;a), i = x,y,z 
and a is the lattice constant. Cyclic permutations yield 
h xz , X z(k) and h yZiyz (k). The U denote effective hopping 
integrals representing the V-O-V hybridization, where 
£o,2j £i,3j and £4 specify the interaction between first, 
second and third neighbors, respectively. For symme- 
try reasons off-diagonal elements arise only between sec- 
ond and third nearest neighbors and are of the form 



-t[Sy 



where Sj 



2 sin(fcja), 



i.e., they vanish at high-symmetry points. Since the 
t[ 2 are very small we neglect the off-diagonal elements 
so that the energy bands are given by = hi y i{k) 1 

with i = xy, xz, yz. The parameters and to ■ ■ ■ £4 can 
be found by fitting the LDA energies at high-symmetry 
points of the Brillouin Zone. 

Fig. 1 (a) shows the t% g bulk bands of SrV03. Accord- 
ing to the predominantly planar electronic structure, the 
bulk density of states Pb{^) exhibits the characteristic 
asymmetric peak related to the van Hove singularity at 
the X point (see Fig. 1(b)). The overall shape of the bulk 
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FIG. 1: (a) Tight-binding fit to LDA t 2g bulk bands of 
SrV03 (Ef = 0). (b) Solid curve: isotropic bulk density of 
states; dashed and dotted curves: local density of d xz , yz and 
d xy states in the first layer of SrVC>3. 



density of states agrees with LDA calculations At 
the surface, the t2 g degeneracy is lifted since only the d xy 
band exhibits strong dispersion within the surface plane 
(the z axis defines the surface normal). To evaluate the 
surface density of states we use a Green's function for- 
malism fl9(l for semi-infinite tight-binding systems. As 
electronic structure calculations for SrVC>3 surfaces are 
not available, we assume the tight-binding parameters to 
coincide with those in the bulk, except for a small surface 
potential to ensure charge neutrality. 

Fig. 1 (b) shows the local density of states p s (w) of the 
d xz ,yz bands for the surface layer of SrVC>3. Its spectral 
weight is reduced at low and high energies but enhanced 
at intermediate energies. Thus, the effective width of 
p s (w) is smaller than that of Pb{<+>) although their total 
widths are identical. In contrast to the d xz ^ yz density, the 
density of d xy states nearly coincides with the bulk den- 
sity, reflecting the planar nature of the ti g states. In the 
deeper layers the local density oid XZtVZ states approaches 
Pb(^) rather quickly, the main effect consisting in an 
oscillatory distortion of the spectral shape. This rapid 
convergence is a consequence of the tight-binding char- 
acter of the t2 9 bands. It would be desirable to carry out 
self-consistent electronic structure calculations for SrVC>3 
surfaces since they should provide a more accurate den- 
sity of states. Nevertheless, the key effect discussed here, 
namely, the preferential narrowing of the local density of 



d X z,yz states, should hold quite generally. 

To interpret the experimental photoemission data we 
evaluate the quasi-particle spectra by accounting for lo- 
cal Coulomb correlations. According to the one-electron 
properties discussed above we are dealing with a non- 
isotropic system where two narrow bands interact with a 
wider band. In the first layer this difference in effective 
band width is largest and it diminishes rapidly towards 
the interior of the system. This situation is reminiscent of 
the layer perovskite Sr2Ru04, which essentially consists 
of Ru sheets containing two narrow d xz . yz bands inter- 
acting with a wide intra-planar d xy band. The peculiar 
feature of this system is the fact the on-site Coulomb 
energy lies between the single-particle widths of the non- 
degenerate t 2g bands: W xz , yz < U < W xy p0| |. In the 
present case the difference between d xz ,yz an d d xy bands 
at the surface is less pronounced so that Wi < U for 
all three bands. Nevertheless, since in the first layer the 
effective width of d xz . yz states is reduced, the influence 
of Coulomb correlations on the surface bands should be 
stronger than on the wider bulk bands. 

The key quantity characterizing the effect of Coulomb 
correlations is the complex self-energy. Since we neglect 
the weak hybridization between ti g bands, S(oj) is di- 
agonal in orbital space. Despite this simplification, a 
full description of correlations near the surface would be 
exceedingly complicated since, in principle, it would re- 
quire a mixed momentum/real space approach in order 
to handle the loss of translational symmetry normal to 
the surface. This could be accomplished using a cluster 
formalism where the semi-infinite system is represented 
via a slab of finite thickness. Unfortunately, the pla- 
nar character of tig orbitals leads to a slowly conver- 
gent local density of d xz%yz states, with many spurious 
van Hove singularities stemming from the quasi- 
one-dimensional hopping parallel to the surface. On the 
other hand, a cluster generalization of the DMFT is fea- 
sible today only for very small cluster size. We therefore 
ignore the momentum variation of the self-energy and as- 
sume that, for a given layer, it depends only on the local 
density of states within that layer ||. The reduced di- 
mensionality at the surface enters the many-body calcu- 
lation via the layer-dependent density of states. In view 
of the local Coulomb interaction this assumption seems 
justified for qualitative purposes. 

Before discussing the quasi-particle spectra derived 
within the DMFT it is useful to compare the frequency 
dependence of the self-energy in the bulk and at the sur- 
face. For convenience we give here the results obtained 
within non-self-consistent second-order perturbation the- 
ory since this provides a clear identification of the vari- 
ous orbital contributions to the diagonal elements 
pfj| . Fig. 2 illustrates the second-order self-energy for 
U = 4.3 cV and J = 0.7 cV. Although the total widths 
of the bulk and surface densities shown in Fig. 1 are the 
same, the effective band narrowing at the surface yields 
a correspondingly narrower distribution of IvaT, xz , yz {uj) 
compared to ImE(,(w) for the isotropic bulk. Because 
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FIG. 2: Second-order self-energy of SrVOa. Solid curves: 
bulk tig self-energy; dashed (dotted) curves: dxz,yz [dxy) self- 
energy for surface layer. For clarity the real part is shifted 
down by 3 eV. The straight line gives u> — ek for eu = —0.5 eV. 



of the orbital mixing via the on-site Coulomb inter- 
action ImS I9 (w) lies between these two distributions. 
From the Kramers-Kronig relations it then follows that 
Re ^xz,yz(u) exhibits a more pronounced minimum be- 
low the band bottom than Re£&(a;). Thus, as indicated 
by the line w — for a typical band energy ej., the surface 
satellite appears at a lower value of U than in the bulk. 
Analogous arguments hold for the upper Hubbard peak. 
This qualitative picture is fully supported by the more 
accurate DMFT results discussed below. 

In the DMFT the self-energy elements Si (us) are func- 
tionals of the bath Green's functions Q^ 1 = G^ + Si. 
The quasi-particle densities of states are related to the 
local Gi via Ni(u>) = —ImGiioS)/^. On-site Coulomb 
correlations are treated using the self-consistent multi- 
band Quantum Monte Carlo method (for a review, see 
Ref. fft). The temperature of the simulation is 12.5 meV 
with 64 imaginary time slices and several runs using 10 5 
Monte Carlo sweeps. The quasi-particle density of states 
is obtained via maximum entropy reconstruction pl]| . 

Fig. 3 (a) shows the bulk quasi-particle density of states 
of SrVC-3 for two Coulomb energies: U — 4.0 eV and 
4.3 eV. The exchange energy is J = 0.7 eV §|. These 
results show that in the bulk U must be larger than 4 eV 
to obtain the satellite observed in photoemission spectra 
jl], ||]. The peak near 2 eV above Ep agrees with in- 
verse photoemission data Although for U = 4eV 
there is considerable correlation-induced band narrowing 
and an emerging satellite shoulder, the increase of U to 
4.3 eV yields an even narrower coherent feature near Ep, 
with the missing weight shifted to the lower and upper 
Hubbard bands. Note that Nb(Ep) — pb(Ep) which fol- 
lows from the local approximation implicit in the DMFT 
for isotropic systems p3| . The most recent photoemis- 
sion data H confirm this result. The spectra shown in 
Fig. 3 (a) qualitatively agree with DMFT-LDA results 
for the 3d 1 perovskite La^Sri-j; Ti0 3 H which exhibits 
a similar ti g bulk density of states. 

The surface quasi-particle spectra for SrVC>3 are shown 
in Fig. 3(b) and (c). The lower Hubbard peak of the 
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FIG. 3: Quasi-particle density of states of SrVOs derived 
from DMFT. (a) bulk t2 g states, (b) surface d xz , y z states, 
(c) surface d xy states. Solid curves: U = 4.3 eV; dot-dashed 
curves: U = 4.0 eV; dashed curves: bare densities. 



dxz.yz states is clearly visible already for U — 4eV be- 
cause of the narrower local density of states in the first 
layer. A larger U shifts the satellite to higher binding 
energies. The comparison with the spectra shown in 
Fig. 3 (a) demonstrates that correlation effects for a fixed 
value of U are stronger at the surface than in the bulk: 
The coherent peak near Ep is narrower at the surface 
and the incoherent satellite feature is more pronounced 
than in the bulk, in agreement with experiment jl], ^, [s| . 

The surface quasi-particle density of d xy states in 
Fig. 3(c) is intermediate between Nb(uj) and N s (ui) for 
d X z,yz- Although there is little single-electron hybridiza- 
tion between ti g bands, the local Coulomb interaction 
mixes them so that the d xy surface spectrum involves 
contributions arising from the more strongly correlated 
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dxz.yz states. Another consequence of the anisotropic sur- 
face self-energy is the correlation-induced charge transfer 
between subbands 0. Here, we find that 0.06 electrons 
are shifted from the d xy to the d xz ^ yz bands. Also, the 
quasi-particle partial densities of states at Ep do not 
need to coincide with the bare partial densities. The cou- 
pling between narrow and wide bands is a genuine multi- 
band effect and underlines the fact that single-particle 
bands in the presence of local Coulomb interactions can- 
not be considered independently. This issue will be con- 
sidered in more detail elsewhere j2f|. 

Note that the many-body reduction of the quasi- 
particle band width is much larger than the surface- 
induced one-electron band narrowing. On the other 
hand, since with U /W » 2 the on-site Coulomb energy is 
not far from the critical value for a metal-insulator tran- 
sition , the band narrowing substantially enhances the 
influence of correlations at the surface. 

It would be interesting to perform angle-resolved pho- 
toemission measurements to determine the correlation- 
induced band narrowing of the SrV03 bands. For 
instance the energy at f should be only a few tenths 
of an eV below Ep instead of 1 eV as predicted by the 
LDA. Accordingly, the true t^g bands should be consid- 
erably flatter than the single-particle bands. Also, mea- 
surements using polarized light could help to separate 
correlation effects in the d xy and d xz ^ yz bands. 

We finally comment on the difference between CaV03 
and SrV03. According to LDA band structure calcu- 
lations for CaRu03 and SrRu03 |25|] , orthorombic dis- 



tortions cause a slight narrowing of the ti g bands and a 
strong broadening of the van Hove singularity above Ep. 
Assuming the same on-site Coulomb energy for the V 4+ 
ions in CaV03 and SrV03, a smaller bulk band width 
of CaV03 implies a transfer of spectral weight from the 
coherent to the incoherent peak. At present, the avail- 
able data g U [fl], |H do not yet fully agree on the 
magnitude of this effect. At the surface, such a transfer 
of spectral weight should be more pronounced because of 
the band narrowing discussed in the present paper. This 
effect is indeed observed in the photoemission spectra 
|^] . For a more detailed analysis it would be important 
to carry out electronic structure calculations for CaVC>3 
in order to study orthorombic distortions and possible 
reconstructions at the surface. 

In summary, we have performed DMFT quasi-particle 
calculations for SrV03 in the bulk and at the surface. 
As a result of the planar nature of the ti g states and the 
surface narrowing of the local density of states on-site 
correlation effects are more pronounced at the surface 
than in the bulk, in agreement with photoemission data. 
The two-dimensional character of d states near Ep is one 
of the hallmarks of transition metal oxides. The surface- 
induced enhancement of correlation effects discussed in 
the present work should therefore be a phenomenon ob- 
servable in many materials. 

I thank A. Bringer for discussions and A.I. Lichtcnstein 
for the QMC-DMFT code. I also thank S. Biermann for 
comments and for bringing Ref. 2 to my attention. 
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